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Annual Summary 

The work performed during the first year of the LOKI Wind Correction Computer 
development program has fallen into four phases, described in detail in the North 
American Instruments, Inc. Quarterly Report series. 

In the first phase, the effect of an arbitrary wind distribution on the LOKI trajec- 
tory during boost was determine by. an impulse response method, and a function was 
derived relating burnout ieee of flight path with the distributed wind along the 
boost trajectory. Complete details of the derivation are included in the First Quart- 
erly Report, Ref. 5. To verify the theoretical influence function, arrangements have 
been made for an spect ieatal sctiparteus at White Sands Proving Ground. The 


Signal Corps is erecting ten poles instrumented with s<nemometers at the launching 


site of the Small Missiles Range, for this purpose. This system is expected to be in 


operation within the next few months. 

In the second phase, the wind influence function was then applied to arbitrary and 
observed wind distributions with a view to establishing the basis for a practical field 
wind measuring system to transmi: iniormaticu to a computer which could perform 


an aiming correction. The studies indicated the feas‘bility of a correction method 
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based on a linear relationship between burnout deviation and wind, suitably measured, 
at a limited nuu.ver of fixed stations. This work is described in the Second Quarterly 
Report. The conclusions are supported by further studies included in this report. 

To supplement immediately the wind data obtained from the literature, which ex- 
tend up to 250 feet altitude (50% of the burnout deviation effect), with measurements 
extending toi, COC feet of altitude (94% of the burnout deviation effect), a balloon-borne 
het wire udeluumeter system was constructed and put into field operation, data being 
gathered inthe Mojave Desertat E] Mirage Airport, and at White Sands Proving Ground 
ai the S::7'! Missiles Range. The equipment is described in the Third Quarterly Re- 
port. A more complete program involving space and time variations of wind within a 
volume significant in the general rocket aiming problem has been set up under the 
spousorship of the Signal Corps. The system invclves the arecten of a number of high 
towers nd the operation of a complex data recording system. It is expected that data 
will become available from this project within the next year. 

At the outset of the computer development program the possible need was recognized 
for an anemometer capabie of fast response, to provide orthogonal wind components 
in terms of electrical outputs. No such instrument was available and a design was 


undertaker based on the drag iorce exerted by the wind stream on a transverse cylinder. 
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The develo, neat of this 2-ccmponent anemometer constituted the third phase and extended 
through most of the year and resulted in a prototype unit suitable for field operation. De- 
tails on the aesign and wind tunnel testing of the anemometer are contained in the Third 
Quarterly Report. 

The final phase of the past year's work, described in this report, has been the design 
ofa complete wind correction computer system based on the influence functionand the re- 


sults of the wind studies. 
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WIND CORRECTION COMPUTER 


Introduction 

The purpose of the computer which is the subject of this report is to correct the aim 
of a rocket launcher (in both azimuth and elevation) for the wind forces experienced by 
the rocket during its burning period. Asaresultof a statictical study of winds within a 
few hundred feet of the ground, it is believed that the effect of these so-called surface 
winds upon the rocket. trajectory can be expressed by a linear relation involving the two 
horizontal wind components at or near the rocket launcher. The details of the wind 
structure study are presented in this report and in am earlier one, Ref. 1, and the ene- 
mometer used in both the wind study and ip the wiud correction computer has been de- 
scribed earlier, Ref. 2. 

The angular deviation at the end of burning due to wind normal to the trajectory, a- 
mounts to apprcximately one mil (thousandth radian) per mile-per-hour; it is tot expect- 
ed that rocket launching will be attempted in winds above 40 miles per hour. Current 
experimental rocket dispersions amountto 8 mils (linear standard deviation), and a wind 


correction standard deviation of about 2 mils thus appears acceptable. 
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Operating Principle of Computer 


The wind correction computer corrects the point of aim of the rocket launcher for 
deviations fromits intended trajectory due to wind. The corrections are made in both 
azimuth and elevation; both corrections are inserted as rotations of the shafts of diff- 
erential selsyns interposed in the transmission lines between the main fire control 


director (M33) and the launcher. The wind correction computer is physically separ- 


RM ee 


ate from the main fire control director and may he remote from it. 

The wind data received from the anemometer is in the form of two amplitude-mod- 
ulated 3000 cps signals; the amplitudes are proportional to v2cos 6 and v’ sin O , 
the northerly and westerly components of the square of the wind velocity; the wind is 2 
assumed to lie in the horizontal plane. Here V is the wind velocity anc © is the wind 
vector heading angle measur: . counter-clockwise relative to true north. Because the 
anemometers do not deliver actual velocity components, the computer must convert 
Vcos © and yin © iuto Vcos © and Vsin © , respectively. Secondly, if the 
anemometer is fixed in orientation, (i.e., not slaved in azimuth to the launcher), it 
is necessary tocarry out the vector rotations which resolve the northerly and westerly 
wind components into down-rarge and cross-range winds; this is achieved in the com- 


puter by sine and cosine potentiometers and operational adders. If A is the azimuth 
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angle of the launcher, the desired quantities are: 
Veos(@ -A) = (Vcos © ) cosA¢ (Vain © ) sin A (1) 
Vain(@ -A) = -(Vcos © ) sinA + (Vsin 6 ) cos A (2) 
We may replace (@ -A) by G, the wind heading relative to the launcher azimuth. 
Since only those wind components normal to the trajectory are effective in produc- 
ing deflections of the trajectory, tie angular rotation duc tc wind may be resolved into 
two angular rotations normal to the trajectory. The horizontal and vertical angular 
deviations W, and W., respectively, may be obtained by the geometric relations, ill- 
ustrated in Fig. 1. Because of the smallness of the angular deviation, the deviation 
resolves in the same manner as the wind force which causes it. The incremental 
angles are: 
“Wy, 2 KVain G (3) 
Wy «KV cos 9 sink (4) 
where E is the corrected angle of launcher elevation and K is the influence cocffi- 


client relating wind (in mph) to angular deviatiun (in mils). 


It should be remembered that the deviations are angular and the incremental vec- 
tors, Wy, Wh, and Wy are chords of a unit sphere. Fig, 2 illustrates the spherical 


trigonometric relations utilized ia deriving the wind corrections in azimuth and elevation. 
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FIG. 1 RESOLUTION OF WIND VECTOR OR ANGULAR DEVIATION 


NORTH AMERICAN INSTRUMENTS INC. 
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FIG. 2 RESOLUTION OF ANGULAR DEVIATION DUE TO WIND 
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By the law of sines, one may write 
sin dA = sin W;/cos E' 
has (5) 


The quantity sin W, may be replaced by W, to within one part in 4,000; sin dA 
cannot be replaced by dA because dA may be as great as 560 mils; thus the equation 
for dA may be written: 

sin dA «= W, sec E' an (:)) 
The correction in elevation, dE, may be written by a second application of the law of 
sines as; 


sin (E¢ W,) - sin E'/cos Wp (7) 


We shall want the expression iv a form not involving E, so we substitute 


E »#» E'e@ dE 
(8) 
with the result: 
dE « -W Wh tan E' 
E «- yt ae ee (9) 
2 
In writing Equation (9), we have made use of the approximations: 
cos (GE# W,) £1 0. 3% (10) 
sin (dE#W,) 5 (E+ Wy) 0.1% (11) 
Ww Z -§ 
cos W, ie —_-— 10 "% (12) 
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Equations (3), (4), (6), and (9) provide tie mathematical basis for the computation wf 

windcorrectiona. Equation (4) involves the angle E while all other relations involve E'. 
If we denote 


W'., = KV cos 6 sin E' (13) 


Wi2 
then (1+ W', cot E')dE = -W'i¢ tan E’ (14) 
2 


Because of the small magnitude of W',, neglect of the term W', cot E' causesanerror 


of less than 4% in dE. Hence, to within 1.6 mils, Equation (9) may be replaced by 
Wh 2 
2 


dE = -Wio¢ tan E! (15) 


the corrected elevation, E, does not appear in any of the computer relations and no 
servo shaft positioning device need be provided for it. If desirable, the computer can 
be arranged to avoid this slighterror. Fig..3indicatesa solution utilizing Equation (14). 
Computer Block Diagram 

Fig. 3 presents a block diagram of the wind correction computer. The five hori- 
zontal dashed lines represent mechanical shafts, positioned by servomotors on the 
right; several linear or functional potentiometers are placed on each shaft; indicating 


dials (and selsyn or potentiometer transmitters, if required) are on the left side of 


the diagram. 
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The manner in which the computer executes the mathematical relations may be seen 
fairly readily. The two anemeometer component aignals Ww sin © and w cos O are 
amplified to auitable level and introduced at the top of Fig. 3. The three linear poten- 
dometers and the quadrature superposition unit are used for extracting the square root 
of V" and dividing each of the input signals tu yield W cos @ and Wain © . These aig- 
nals are combined by means of the azimuth sine-coasine potentiometers and aubsequent 
adding amplifiers to produce the components of the correction relative to the instantaneous 
launcher heading, W cos — and W ain @. 

Several functional potentiometers are employed to modulate the V cos G and V sing 
signals by sin E', cos E', ‘al ee E', the outputs from these functional potentlometers 
are combined in several ways by adding amplifiers and the resulting signals are balanced 
by output servos against the signals derived from other potentiometers maunted on the 
output shafts. 

Square-Root Techniques 

The poteutiometer cascade utilizedin the computer schematic of Fig. 3 is rather simple 
in principle, but may lead to significant dynamic error in operation due to the slowness 
of the servo-balancing as compared to wind fluctuations. An analysis of the dynamic 


error due to this cause has not yet been carried out; however, it is apparent that the error 
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x 
may amount to several percent unless the input signals are damped to correspond to the 
servo capabilities. 

Damping may be achieved by means of an orifice inthe anemometer (up to about 15 sec- 
onds time constant) or by filter networks placed in dc. sections of the computer; exper- 
imentally, the latter method is more appealing because of the possibility of awitching 
damping values. However, it is convenient to avoid the troubles of d.c. amplification 
in the computer by using the 3000 cps signals throughout; this favors mechanical damping 
To lessen the importance of this controversy, it is contemplated to make the experimental 
computer convertible insofar as possible from a.c. tod.c. signal operation; this is the 
reasonreactive resoivers are not used. If d.c. operation is employed, a.c. signals are 
none the less required through the quadrature circuits, and the rectification must take 
place after the square root has been extracted. 

To eliminate the possible aia of dynamic error in aquare rooting, two non-servo- 
square-root circuits have been investigated. One type of square-rooter utilizes the non- 
linear properties of Thyrite, while a second type employs suppressor controlled pentodes 
ina feedback circuit to extract the square root; simplified breadboard computers have 
beex built utilizing both approaches. The two methods of square rooting are treated 


further in the following section. 
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Thyrite Type Computer , 


The computer utilizing Thyrite as a square root extracting device combines the 
signals V“ cos @ and V’sin © inphase quadrature to produce a signal of ampli- 
tude V". The Thyrite square root extractor delivers d.c. voltages proportional to 
+1 Vi and - {v| , which are applied to the two potentiometers. The output angular 
positions representing the azimuth and eicvation corrections — positioned by servo- 


mechanisms effecting the following relations: 


&. > mame (ine oe + Oo btn OTR otvpellges i 


lv|sin OA w KV“ Bing sec E' (16) 


| 

| V| siv AE -~ KV" cos @ sin E! (17) 

The computer requires the solution for which Iv| yO; for this reason [v| is : 

constrained to be significanily greater than zero at all times, regardless of wind. : 

The use ofa variable voltage lv across the servo potentiometers leads to a system 

in which the servo loop gain is variable; as is wiil xGowa, this situation is usually 
conducive to instability under some of the desired operating conditions. This situ- 
ation might be tolerated in this mode of operation because of the convenience and 


economy of effecting both the square root extraction and the output positioning with 


the same servo. 


— 
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ite Characteristics 

The non-linear characteristics of Thyrite, a material manufactured by the General 
Electric Company for use in lightning arrestors, are demonstrated by the plot of Fig. 4; 
the graph shows the voltage across a small Thyrite element as a function of the current 
passing through it. The slight departure from the desired square-root dependence is 
almost entirely corrected by the use of a trimming resistor, in this case about 300 ohms. 

A vary simple Thyrite driver circuit was tried, which requires only moderate 
voltages; the circuit is shown in Fig. 5 and its performance is represented by Fig. 7. 
The 5% accuracy achieved in this trial circuit can undoubtedly be improved to corres- 
pond to the 2-3% result of the Thyrite alone. 
Suppressor Co 

To overcome the necessity for d.c. servo amplificationand toavoidthe. hazards of 
variable servo-loop gain, some tests have been made of aquare root extraction by means 
of successive suppressor controlled variable gain amplification stages. The Suppressor 
Type Computer utilizes this type of circuit for square root extraction. The computer 
shown in Fig. 3 is different from the Suppressor Type only in its use of two successive 


potentiometer stages for root extraction rather than suppressor controlled amplification 
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stages. This type of computer, whether it uses suppressor amplifiers or potentiometers 
for root extraction, has the advantage that no d.c. amplification is required; only simple 
phase sensitive rectification is required at the output stages of the servo amplifiers. 
Suppressor Amplifier Root Extraction 
The principle of square root extraction with a feedback circuit incorporating suppressor 
amplifiers is outlined in Fig. 6. The amplifiers labelled "8" are stages employing, for 
example, 6AS6 pentodes in which the gain can be controlled readily by the suppressor grid 


bias. Assuming linear characteristics in the range of operation, the gain of each stage is 


“leVle4K + &B aAx 
sh 2a A5x (18) 


where K is the gain of the suppressor amplifier at zero bias, § is the change of gain per 
volt of bias, aA is the constant portion of the loop gain (subsequent to the suppressor am- 
plifiers and including the rectification), and X is the input signal (RMB volts). A has the 
units d.c. volts per RMB voli. B is the d.c. voltage that is compared with the output of 


the rectifier. The output of the first suppressor stage is, therefore: 


Y » GX x 1 eee ee oe 
w= (HrbBE (1. Bad 6(K+ 5 B)x ) 2aA 6 = 
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Typical values of the parameters are: 
K = 65 
& =6.2 per volt 
9A we §7 
B = 50 volts d.c. 


The calculated behavior of this system shows an error equivalent to less than 0.5 mil; 


oer 6 ee as din Sin she, 


experimental results are not yet available. 

The suppressor amplifier aquare root scheme assumes that the gains of the several 
suppressor amplifiers vary in unison; to the extent that this is required it constitutes a 
strong reasonagainst the use of the scheme in an operation computer. However, it can be 
shown that compensation for differences in tube characteristics can be achieved in large 


measure by interchanging tubes so that, for example 


and adjusting the loop gain so that, approximately 
{ 
1Gpx =o [1+ oe 
1G X C lle (x) (21) | 
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WIND STUDIES FOR LOKI 

In an earlier report of this project, Ref. 1, a series of vertical wind profiles was 
giudied to determine the feasibility of a wind correction computer for LOKI, based on 
a linear relationship between burnout deviation and wind ai a selected level. Only a 
limited number of wind profiles was inyolved in the comparigon. However, significant 
correlations were found between computed burnout deviation and corresponding single 
level windmeasurements. The effects of altitude and smoathing were examined, com - 
parisons being made with wind level and time averaged velocity as parameiers. 

The wind data used was that of Sherlock and Stout for winter storms in Michigan, 
Refs. 3 and 4, andwas obtained with a 250-foot tower instrumented with fest response 
pressure piate anemometers at 25-foot intervals, It was learned from the authors that 
the published information was only a smal} part of a large total which had been reduced 
from the original oscillograph reeords to punched ISM cards. Several thousand of 
these cards were available; duplicates were ohiained and these were processed by the 
Los Angeles IHM Service Bureau to extend the initial investigation, The results of the 
IBM computations which are presented in this report, confirm the initiai findings and 


permit a selection of optimum altitude and time average through the use of a wider range 
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of these parameters. Also comparisons were made using combinations of wind values at 
two altitudes. These appear to offer little improvement over appropriately smoothed 
single altitude winds at levels above 50 feet. 
Discusston 

The function relating LOKI deviation at burnout with the distributed wind along the 
trajectory was derived in Ref, 5 and is shown in Fig. 8 asa plot between wind influence 
coefficient deviation of flight path at burnout from direction at launching per unit cross- 
wind per unit distance increment along trajectory) and distance of rocket from launcher. 

The cumulative effect of a constant crosswind during boost is shown in Fig. 9. The total 
deviation for unit crosswind is .8 mil, 50% being accumuwiated in the first 200 feet of tra- 
jectory and 95% in the first 1,000 feet. 

The results of the preliminary investigation are tabulated in Fig, 10from Ref. 1. The 
burnout deviations for 150 cases were calculated by integrating the wind observations 
weighted with the influence function and apply to a theoretical vertical shot. Since the 
wind observations extend to 250 feet, the deviation is about 50% of the total that would be 
obtained if the integration were carried to burnout. 

The work discussed inthis report extends the number of cases from 150 to 4, 000 using 


the same altitudes, 25, 50, 100, and 150 feet with the addition of the mean of the 25 and 
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Wind 4-Second Average Wind | 
Wind Deviation at Burnout Deviation at Burnout i 


FIG, 10 SUMMARY OF RESULTS OF PRELIMINARY WIND STUDIES FOR LOKI 
(Ref, 1) 
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100-foot winds, and the mean of the 50 and 150-foot winds. Also the 1 and 4-second 
wind averages in the initial survey are augmented by the 2 and 8-second averages. 

The basic wind data as received from the University of Michigan, was in the form 
of punched IBM cards and apply to the storms of Apr, 28, 1931 and Jan. 19, 1933. 

All records were originally taken continuously on oscillograph paper, information 
being transferred from these records onto the cards. The card code identification 
is shown in Fig. 11. 

The first vertical column at the left identifies the run number. Each storm sample 
was divided into runs of various lengths, the oscillograph being operated intermittent- 
ly to record during periods of major gust activity. Many runs were found to be bro- 
ken, due to missing card data, and in order to provide continuous time sequences 
for analysis, these runs were broken into sub-groups and labelled episodes. 

The second coluin identifies the length of time interval over which the velocity 
observations were averaged. Averages from 1/4 sec. to 10 sec. were available; 

however, the averages used in all computations were based on the 1-second cards. 

The third column shows the number of seconds after the start of a run, and col- 
umns § to 16 show the wind speeds in miles per hour at each of the observation sta- 
tions. The ten stations on the tower were at 25-foot intervals and are {dentified by 


numbers 1 to 10 reading from the bottom. 
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a) 6= 
Statistical comparisons were made for each run separately to determine the correlation, 
coefficient, the slope of the burnout deviation vs. wind speed regression line and the 90% 
confidence limits for the burnout deviation. In calculating the burnout deviation sequences, 
the wind values at each altitude were assumed to describe the wind extendirg from the acja- 
cent lower altitude using weighting factors determined from the influence function. The 


results as tabulated by IBM are presented in the Appendix. 


The formulae used in the computations follow: 


Correlation coefficient, r a N& dA-2A A-= AS d z 
[sat ceA)*] [N= 5*-(f 37} : 


2 a 2 Ve 
90% Confidence Limits, Sygq = LGS [1-at TINE §°- CE $) J 


MN 


NEEA- E5EA 
NE§'‘=- (= &)* 


Regression Line Slope, b = 


(Regression Line passes through A > ) 
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Burnout Deviation, é vk 2.4 Mijn , mils 


where : 
K ie the episode number } 
j is the time in seconds after the astact of a run a 
m is the index for moving averages, m:=1l, 2, 4, 8 
W ia the velocity in miles per hour 
A jk are the wind weighting factora (see table below) 
K is a conversion constant, K = 88/60 
i is the height af observation given in the following table: | 
: 
i He Wind w factor | 
1 25 . 050 
2 50 . 050 
3 75 . 930 
4 100 . 035 
5 125 . 030 
6 150 046 | 
7 175 , 050 
8 200 . 640 
9 250 . 066 
p is a specified value or combinution of i's as given in the following table: 
-: Ww W, 
“Py ty %% z= 46 
W W 
i ¥ i 
~, * & Mp Se ee 
WwW 
2 
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Results are summarized in Figs, 12 through 14,which show the range and average of 
the three aignificant parameters plotted as vertical bars and each of the combinations 
of anemometer height and time average. Runs 17, 19, 34, and 35 have heen excluded 
in certain cases where missing data have resulted in small sample size and resulting 
rangeextremea. The unweighted arithmetical mearis shown on ench har. Fir.J2 shows 
the effect of progressive iucrease of time average on correlation coefficient of burnout 
deviation vs. wind speed for progressively increasing anemometer heights from 25 to 
150 feet and also for time average 25 to 100-foot winds and 50-150-foot winds. Similar 
progressive sequences are shown for the regression line slope in Fig,13 and 90% con- 
fidence limits in Fig. 14. 

It ig evident that the 2-second average offers improved correlation and reduced 
spread among runs under the 1 and 4-second averages. Further, a pronounced degrada- 
tion is observed with the 8-second average. 

In progressing upward in altitude a marked improvement occurs in going from 50 to 
100 feet with subsequent slight improvement at 150 feet. Best resulis are observed in 
the 50-150 foot mean wind sequences. However, the {mprovement over the 100-foot 
wind would hardly warrant the use of 2 anemometers ina wind correction computing 


system. 
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(University of Michigan Wind Data) 
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FIG, 13 EFFECT OF ALTITUDE AND AVERAGE ON SLOPE OF REGRESSION LINE, 
LOKI BURNOUT DEVIATION VS, WIND SPEED 
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The present results are in substantial agreement with those obtained in the prelim-~ 
inary analysis of Ref. 1, amd indicate the best choice of altitude and smoothing that 
wouldapply to vertical trajectories under the condition that would obtain at the Michigan 
test site in winter storm conditions. 

For ananemometer height of 100 feet and 2~second average, the correction for wind 
effect in the first 250 feet of a vertical trajectory (about 50% of the total) could be ex- 


pressed by the relation: 


0.48W + 4 


-~&ec 


where § = burnout deviation correction, mils 


w! 100 foot 2-second average wind velocity, mph, 


$ 
ae atria het 


| 
. CONFIDENTIAL 
5 
the 90% confidence limits would be + 2 mils. 
Supplementary data obtained from hot wire surveys up to 1000 feet at White Sands 
Proving Ground and El Mirage Airport, using equipment described in Ref. 3, is now 
being analyzed and should permit the establishment of a tentative wind correction cor- 
rection equationfor the entire boost phase. i 
The effects of horizontal separation between the anemometer and launcher, elevation : 
angle of launching, terrain and weather conditions, must also be considered in the de- 
sign of a tactical system. These will be investigated in the Signal Corps sponsored 
project which is now being condu::ted by North American Instruments, Inc. 
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APPENDIX 
STATISTICAL RESULTS OF LOKI BURNOUT DEVIATION 
VS. WIND VELOCITY COMPARISON 

The tables which follow present the results of the statistical computations described 
above, p. ~5, and are direct photographs of the original IBM records. The accuracy is 
no better than 3 figures, the additional figures were allowed to appear as a matter of 
convenience in machine operation. 


The identification code for the various columns as given on p. 15 are repeated below: 


p = height of observation code number 
m - index for moving average, seconds 
k w» episode number 


N = number of observations 

Ag # mean burnout deviation, mils, for episode 
A =» mean wind velocity, mph, for episode 

b = regression line slope, mils/mph 

By 90 z= 90% confidence limits, mils 


r » correlation coefficient 
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